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I review recent analytical and numerical advances in the study of non-equilibrium quark- 
gluon plasma physics. I concentrate on studies of the dynamics of plasmas which are locally 
anisotropic in momentum space. In contrast to locally isotropic plasmas such anisotropic plas- 
! mas have a spectrum of soft unstable modes which are characterized by exponential growth of 
transverse (chromo)-magnetic fields at short times. Parametrically the instabilities provide the 
fastest method for generation of soft background fields and dominate the short-time dynamics 
■ of the system. 
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! 1. Introduction 

Ph! One of the most important open questions emerging from the Relativistic Heavy-Ion Col- 

lider (RHIC) program is how quickly, and by what process, the stress-energy tensor of the 
high-density QCD matter produced in the central rapidity region becomes isotropic. This is 
important, for example, to determine the maximal temperature achieved in such collisions, 
which will also be pursued in the near future at even higher energies at the CERN Large 
Hadron Collider (LHC). One of the chief obstacles to thermalization in ultrarelativistic heavy- 
ion collisions is the rapid longitudinal expansion of the central rapidity region. If the matter 
expands too quickly then there will not be sufficient time for its constituents to interact and 
thermalize. In the absence of interactions, the longitudinal expansion causes the system to be- 
come much colder in the longitudinal direction than in the transverse directions, corresponding 
to (p 2 L ) <C {pj) in the local rest frame. One can then ask how long it would take for interactions 
to restore isotropy in momentum space. 

High-energy heavy-ion collisions release a large amount of partons from the wavefunctions of 
the colliding nuclei. Partons with very large transverse momenta originate from high-Q 2 hard 
interactions, while partons with transverse momenta below the so-called saturation momentum 
Q s (given by the square root of the color charge density per unit area in the incoming nuclei) are 
much more abundant and are better viewed as a classical non-Abelian field [ EJ 121 E| . At high 
energies, Q s ^> Aqcd sets a semi-hard scale which allows for a weak-coupling treatment of the 
early-stage dynamics. However, due to the large occupation number of the soft modes below Q s 
(the classical color field has strength ~ l/flO; the problem is nevertheless non-perturbative. 
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The evolution following the initial impact was among the questions which the so-called 
"bottom-up scenario" [Hj attempted to answer. For the first time, it addressed the dynamics 
of soft modes ( "fields" ) with momenta much below Q s coupled to the hard modes ( "particles" ) 
with momenta on the order of Q s and above. However, it has emerged recently that one of the 
assumptions made in this model is not correct. The debate centers around the very first stage of 
the bottom-up scenario in which it was assumed that (a) collisions between the high-momentum 
(or hard) modes were the driving force behind isotropization and that (b) the low-momentum 
(or soft) fields act only to screen the electric interaction. In doing so, the bottom-up scenario 
implicitly assumed that the underlying soft gauge modes behaved the same in an anisotropic 
plasma as in an isotropic one. However, it turns out that in anisotropic plasmas, within the so- 
called "hard loop approximation" (see below), the most important collective mode corresponds 
to an instability to transverse magnetic field fluctuations [ U\ • Recent works have shown 
that the presence of these instabilities is generic for distributions which possess a momentum- 
space anisotropy [ |H1 El HHj and have obtained the full hard-loop (HL) action in the presence 
of an anisotropy [lllj. Another important development has been the demonstration that such 
instabilities exist in numerical solutions to classical Yang-Mills fields in an expanding geometry 

Recently there have been significant advances in the understanding of non-Abelian soft-field 
dynamics in anisotropic plasmas within the HL framework [El El El- The HL framework is 
equivalent to the collisionless Vlasov theory of eikonalized hard particles, i.e. the particle tra- 
jectories are assumed to be unaffected (up to small-angle scatterings with 9 ~ g) by the induced 
background field. It is strictly applicable only when there is a large scale separation between 
the soft and hard momentum scales. Even with these simplifying assumptions, HL dynamics for 
self-interacting gauge fields is non-trivial due to the presence of non-linear interactions which 
can act to regulate unstable growth. These non-linear interactions become important when 
the vector potential amplitude is on the order of A non _ Abelian ~ Ps/g ~ (9Ph)/g, where ph is 
the characteristic momentum of the hard particles and p s is the characteristic soft momentum 
of the fields. In QED there is no such complication and the fields grow exponentially until 
^Abeiian ~ Ph/ g a ^ which point the hard particles undergo large-angle deflections by the soft 
background field invalidating the assumptions underpinning the hard-loop effective action. 

Recent numerical studies of HL gauge dynamics for SU(2) gauge theory indicate that for 
moderate anisotropies the gauge field dynamics changes from exponential field growth indica- 
tive of a conventional Abelian plasma instability to linear growth when the vector potential 
amplitude reaches the non-Abelian scale, A non _Abciian ~ Ph [Ellin]- This linear growth regime 
is characterized by a cascade of the energy pumped into the soft modes by the instability to 
higher- momentum plasmon-like modes [El El- These results indicate that there is a funda- 
mental difference between Abelian and non-Abelian plasma instabilities. However, even with 
this new understanding the HL framework relies on the existence of a large separation between 
the hard and soft momentum scales by design. One would like to know what happens when 
the scale separation between hard and soft modes is not very large or when the initial fields 
have large (non-linear) amplitudes which is seemingly the situation faced in real experiments. 
In this case one is naturally led to consider instead numerically solving the Wong- Yang-Mills 
(WYM) equations [EDI HQ ■ 
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2. Anisotropic Gluon Polarization Tensor 

In this section I consider a quark-gluon plasma with a parton distribution function which is 
decomposed as 

f(p) = 2Nf (n(p) + n(p)) + AN c n g (p) , (1) 

where n, n, and n g are the distribution functions of quarks, anti-quarks, and gluons, respectively, 
and the numerical coefficients collect all appropriate symmetry factors. Using the result of 
Ref. [ Ej the spacelike components of the high-temperature gluon self-energy for gluons with 
soft momentum (k ~ gT) can be written as 

where the parton distribution function /(p) is, in principle, completely arbitrary. In what 
follows we will assume that /(p) can be obtained from an isotropic distribution function by the 
rescaling of only one direction in momentum space. 

In practice this means that, given any isotropic parton distribution function f\ so (p), we can 
construct an anisotropic version by changing the argument of the isotropic distribution function, 

/(p) = VI + £ /iso (\/p 2 + £(p • n) 2 )) where the factor of y/1 + £ is a normalization constant 
which ensures that the same parton density is achieved regardless of the anisotropy introduced, 
n is the direction of the anisotropy, and £ > —1 is an adjustable anisotropy parameter with 
£ = corresponding to the isotropic case. Here we will concentrate on £ > which corresponds 
to a contraction of the distribution along the n direction since this is the configuration rele- 
vant for heavy-ion collisions at early times, namely two hot transverse directions and one cold 
longitudinal direction. 

Making a change of variables in (J2J) it is possible to integrate out the |p|-dependence giving [ 

1 

WHulk 9)-u 2 f —v* v l + Zt^L (gSi + v ' kl ) ( 3 ) 

where cos8 n = k n and /j, 2 = \/l + £ m 2 D > 0. The isotropic Debye mass, mp, depends on / iso . 
In the case of pure-gauge QCD with an equilibrium / iso we have mp = gT. 

The next task is to construct a tensor basis for the spacelike components of the gluon self- 
energy and propagator. We therefore need a basis for symmetric 3-tensors which depend on 
a fixed anisotropy 3- vector n % with n 2 = 1. This can be achieved with the following four 
component tensor basis: A l i = 8^ — k % W /k 2 , B 1 ^ = k l k 3 /k 2 , C % i = ^h^/h 2 , and D % i = 
k % hi + Wh % with h l = A^n 3 . Using this basis we can decompose the self-energy into four 
structure functions a, ft, 7, and 5a.sIl = aA + j3~B + 'yC + 5~D. Analytic integral expressions 
for a, j3, 7, and 5 can be found in Ref. [|Hj and [ITUj. 

3. Collective Modes 

As shown in Ref. [ |H] this tensor basis allows us to express the propagator in terms of the 
following three functions 

A-\K) = k 2 -uj 2 + a, 
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In Figure Q I plot the instability 
growth rates, r a and T_, clS db func- 
tion of wave number for £ = 10 and 
6 n = tt/8. Note that both growth 
rates vanish at k = and have a 
maximum T* ~ yt//10 at A;* ~ /x/3. 
The fact that they have a maximum 
means that at early times the sys- 
tem will be dominated by unstable 
modes with spatial frequency 
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Figure 1. Instability growth rates as a function of wave number for £ = 10 and 6 n = ir/i. 



A ± \K) = cu 2 -nl, 

where 2Q 2 ± = Cl 2 ± ^O 4 - 4((a + 7 + fc 2 )/? - Pn 2 5 2 ) and H 2 = a + (3 + 7 + fc 2 . 

Taking the static limit of these three propagators we find that there are three mass scales: 
m± and m a . In the isotropic limit, £ — > 0, m 2 = m 2 = and m 2 ^ = m 2 D . However, for £ > 
we find that m 2 < for all \9 n \ and m 2 < for all \9 n \ < it /A. Note also that for £ > 

both m 2 and m 2 have there largest negative values at 9 n = where they are equal. 

The fact that for £ > both m 2 and m 2 _ can be negative indicates that the system is unstable 
to both magnetic and electric fluctuations with the fastest growing modes focused along the 
beamline (6 n = 0). In fact it can be shown that there are two purely imaginary solutions to each 
of the dispersions relations A~ l (K) = and Al 1 (ii') = with the solutions in the upper half 
plane corresponding to unstable modes. We can determine the growth rate for these unstable 
modes by taking w — > iT and then solving the resulting dispersion relations for T(k). 

4. Discretized Hard-Loop Dynamics 

It is possible to go beyond an analysis of gluon polarization tensor to a full effective field 
theory for the soft modes and then solve this numerically The effective field theory for the soft 
modes that is generated by integrating out the hard plasma modes at one-loop order and in the 
approximation that the amplitudes of the soft gauge fields obey |p|/<7 is that of gauge- 

covariant collisionless Boltzmann-Vlasov equations [1221- in equilibrium, the corresponding 
(nonlocal) effective action is the so-called hard-thermal-loop effective action which has a simple 
generalization to plasmas with anisotropic momentum distributions [EJ. For the general non- 
equilibrium situation the resulting equations of motion are 

F^(A)v» = [vD{A)]W„{x;v), (4) 

where / is a weighted sum of the quark and gluon distribution functions [ 11 and v M = p M /|p| = 
(l,v). 
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Figure 2. Comparison of the energy transferred from hard to soft scales, |£(HL)|, for 3+1- 
dimensional simulations with N w = 20, 100, 200 on 96 3 , 88 3 , 69 3 lattices. Inset shows late-time 
behavior on a linear scale. 

These equations include all hard-loop resummed propagators and vertices and are implicitly 
gauge covariant. At the expense of introducing a continuous set of auxiliary fields Wp(x\ v) the 
effective field equations are also local. These equations of motion are then discretized in space- 
time and v, and solved numerically. The discretization in v-space corresponds to including 
only a finite set of the auxiliary fields Wp(x; v$) with 1 < % < N w . For details on the precise 
discretizations used see Refs. [ IT5| ITo]. 

4.1. Discussion of Hard-Loop Results 

During the process of instability growth the soft gauge fields get the energy for their growth 
from the hard particles. In an abelian plasma this energy grows exponentially until the energy 
in the soft field is of the same order of magnitude as the energy remaining in the hard particles. 
As mentioned above in a non-abelian plasma one must rely on numerical simulations due to the 
presence of strong gauge field self-interactions. In Fig. |2]I have plotted the time dependence of 
the energy extracted from the hard particles obtained in a 3+1 dimensional simulation of an 
anisotropic plasma initialized with very weak random color noise [ To] . As can be seen from 
this figure at m^t ~ 60 there is a change from exponential to linear growth with the late-time 
linear slope decreasing as Nw is increased. 

The first conclusion that can be drawn from this result is that within non-abelian plasmas 
instabilities will be less efficient at isotropizing the plasma than in abelian plasmas. However, 
from a theoretical perspective "saturation" at the soft scale implies that one can still apply 
the hard-loop effective theory self-consistently to understand the behavior of the system at 
late times. Note, however, that the latest simulations [ EE] have only presented results for 
distributions with a finite 0(1— >10) anisotropy and these seem to imply that in this case the 
induced instabilities will not have a significant effect on the hard particles. As a consequence due 
to the continued expansion of the system the anisotropy will increase. It is therefore important 
to understand the behavior of the system for more extreme anisotropies. Additionally, it would 
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be very interesting to study the hard-loop dynamics in an expanding system. Naively, one 
expects this to change the growth from exp(r) to exp( v /r) at short times but there is no 
clear expectation of what will happen in the linear regime. Significant advances in this regard 
have occurred recently [ [T7j . These initial results apply to Abelian plasmas and it would be 
very interesting to incorporate expansion in collisionless Boltzmann-Vlasov transport in the 
hard-loop regime and study the late-time behavior in the non-abelian case. 



5. Wong-Yang-Mills equations 

It is also to possible to go beyond the hard-loop approximation and solve instead the full 
classical transport equations in three dimensions [ I21j . The Vlasov transport equation for hard 
gluons with non- Abelian color charge q a in the collisionless approximation are [ I2H1 E^] , 

p^d, - gq a F^d; - gf ab cAlq c d qa }f(x, P , q) = . (5) 

Here, f(t, x, p, q a ) denotes the single-particle phase space distribution function. 

The Vlasov equation is coupled self-consistently to the Yang- Mills equation for the soft gluon 
fields, 

DfiF ^ = r = g J jSPdqqV"f(t,x,j>,q) , (6) 

with f M = (l,p/p). These equations reproduce the "hard thermal loop" effective action near 
equilibrium [ |23 I2H EI] ■ However, the full classical transport theory ()5I6|) also reproduces 
some higher n-point vertices of the dimensionally reduced effective action for static gluons [ 
I2H] beyond the hard-loop approximation. The back-reaction of the long-wavelength fields on 
the hard particles ("bending" of their trajectories) is, of course, taken into account, which is 
important for understanding particle dynamics in strong fields. 

Eq. (j3J) can be solved numerically by replacing the continuous single-particle distribution 
/(x, p, q) by a large number of test particles: 

/(*, P,q) = ^-E ^ ~ *0O) (2*) 3 <f(P " Pi(t)) 5(q a - qf (t)) , (7) 

^ test ; 

where Xi(t), Pi(t) and cff(t) are the coordinates of an individual test particle and Y test denotes 
the number of test-particles per physical particle. The Ansatz flj|) leads to Wong's equations [ 



~dt 
dpi 

dt 



g q? (E a + v, x B a ) , (9) 



J av = T^E<?>^(*-*i(t))- (ii) 

J* test 4 
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Figure 3. Leftmost figure shows the time evolution of the field energy densities for SU(2) gauge 
group resulting from a highly anisotropic initial particle momentum distribution. Simulation 
parameters are L = 5 fm, phard = 16 GeV, g 2 n g = 10/fm 3 , = 0.1 GeV. Right two panels 
show the corresponding fourier transforms of the electric and magnetic fields at different times 
which are indicated in the legend. 

for the i-th test particle. The time evolution of the Yang-Mills field can be followed by the stan- 
dard Hamiltonian method [|2H] in A = gauge. For details of the numerical implementation 
used see Ref. [I2T]. 

In Fig. El I present the results of the three-dimensional simulations published in Ref. [ |^] . 
The leftmost panel shows the time evolution of the field energy densities for SU(2) gauge group 
resulting from a highly anisotropic initial particle momentum distribution. The right two panels 
show the corresponding fourier transforms of the electric and magnetic fields at different times 
which are indicated in the legend. 

The behavior shown in Fig. |3] indicates that the results obtained from the hard-loop simula- 
tions and direct numerical solution of the Wong- Yang-Mills (WYM) equations are qualitatively 
similar in that both show that for non-abelian gauge theories there is a saturation of the energy 
transferred to the soft modes by the gauge instability. Additionally, as can be seen from the 
fourier transforms in right two panels of Fig.|3]the saturation is accompanied by an "avalanche" 
of energy transferred to soft field modes to higher frequency field modes with saturation oc- 
curring when the hardest lattice modes are filled. A more thorough analytic understanding of 
this ultraviolet avalanche is lacking at this point in time although some advances in this regard 
have been made recently [ 80 . Additionally, since within the numerical solution of the WYM 
equations the ultraviolet modes become populated rapidly this means that the effective theory 
which relies on a separation between hard (particle) and soft (field) scales breaks down. This 
should motivate research into numerical methods which can be used to "shuffle" field modes to 
particles when their momentum becomes too large and vice-versa for hard particles. This is a 
difficult task but work on such algorithms is in progress. Hopefully, using such methods it will 
be possible to simulate the non-equilibrium dynamics of anisotropic plasmas in a self-consistent 
numerical framework which treats particles and fields and the transmutation among these two 
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types of degrees of freedom smoothly. 
6. Outlook 

An obviously important question is whether quark-gluon plasma instabilities and/or the 
physics of anisotropic plasmas in general play an important phenomenological role at RHIC or 
LHC energies. In this regard the recent papers of Refs. [ EH E21 EH1 El] provide theoretical 
frameworks which can be used to calculate the impact of anisotropic momentum-space distri- 
butions on observables such as jet shapes and the rapidity dependence of medium-produced 
photons. 

REFERENCES 

1. A.H. Mueller, Nucl. Phys. A715 (2003) 20, hep-ph/0208278. 

2. E. Iancu and R. Venugopalan, (2003), hep-ph/0303204. 

3. L. McLerran, Nucl. Phys. A752 (2005) 355. 

4. R. Baier et al., Phys. Lett. B502 (2001) 51, hep-ph/0009237. 

5. S. Mrowczynski, Phys. Lett. B314 (1993) 118. 

6. S. Mrowczynski, Phys. Rev. C49 (1994) 2191. 

7. S. Mrowczynski, Phys. Lett. B393 (1997) 26, hep-ph/9606442. 

8. P. Romatschke and M. Strickland, Phys. Rev. D68 (2003) 036004, hep-ph/0304092. 

9. P. Arnold, J. Lenaghan and G.D. Moore, JHEP 08 (2003) 002, hep-ph/0307325. 

10. P. Romatschke and M. Strickland, Phys. Rev. D70 (2004) 116006, hep-ph/0406188. 

11. S. Mrowczynski, A. Rebhan and M. Strickland, Phys. Rev. D70 (2004) 025004, hep- 
ph/0403256. 

12. P. Romatschke and R. Venugopalan, (2005), hep-ph/05 10292. 

13. P. Romatschke and R. Venugopalan, Phys. Rev. Lett. 96 (2006) 062302, hep-ph/0510121. 

14. P. Romatschke and R. Venugopalan, (2006), hep-ph/0605045. 

15. P. Arnold, G.D. Moore and L.G. Yaffe, Phys. Rev. D72 (2005) 054003, hep-ph/0505212. 

16. A. Rebhan, P. Romatschke and M. Strickland, JHEP 09 (2005) 041, hep-ph/0505261. 

17. P. Romatschke and A. Rebhan, (2006), hep-ph/0605064. 

18. P. Arnold and G.D. Moore, Phys. Rev. D73 (2006) 025006, hep-ph/0509206. 

19. P. Arnold and G.D. Moore, Phys. Rev. D73 (2006) 025013, hep-ph/0509226. 

20. A. Dumitru and Y. Nara, Phys. Lett. B621 (2005) 89, hep-ph/0503121. 

21. A. Dumitru, Y. Nara and M. Strickland, (2006), hep-ph/0604149. 

22. J.P Blaizot and E. Iancu, Phys. Rept. 359 (2002) 355, hep-ph/0101103. 

23. S.K. Wong, Nuovo Cim. A65 (1970) 689. 

24. U.W. Heinz, Phys. Rev. Lett. 51 (1983) 351. 

25. PF. Kelly et al., Phys. Rev. Lett. 72 (1994) 3461, hep-ph/9403403. 

26. PF. Kelly et al., Phys. Rev. D50 (1994) 4209, hep-ph/9406285. 

27. J.P. Blaizot and E. Iancu, Nucl. Phys. B557 (1999) 183, hep-ph/9903389. 

28. M. Laine and C. Manuel, Phys. Rev. D65 (2002) 077902, hep-ph/0111113. 

29. J. Ambj0rn et al., Nucl. Phys. B353 (1991) 346. 

30. A.H. Mueller, A.I. Shoshi and S.M.H. Wong, (2006), hep-ph/0607136. 

31. P. Romatschke and M. Strickland, Phys. Rev. D69 (2004) 065005, hep-ph/0309093. 



9 



32. P. Romatschke and M. Strickland, Phys. Rev. D71 (2005) 125008, hep-ph/0408275. 

33. B. Schenke and M. Strickland, (2006), hep-ph/0606160. 

34. P. Romatschke, (2006), hep-ph/0607327. 



